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Abstract 

Vector spherical wavefunctions were derived in closed-form to represent time-harmonic electromagnetic 
fields in an orthorhombic dielectric-magnetic material with gyrotropic-like magnetoelectric properties. 
These wavefunctions were used to formulate the T matrix for scattering by a three-dimensional object 
composed of the chosen material. Furthermore, a closed-form, coordinate-free expression of the dyadic 
Green function for the chosen material was derived. Expressions ascertained for the singularity behavior 
will be useful for formulating volume integral equations for scattering inside the chosen material. A 
bilinear expansion of the dyadic Green function was obtained in terms of the derived vector spherical 
C/2 ■ wavefunctions. 

x, 

1 Introduction 

Closed-form solutions of the vector differential equation 

V x [V x F(r)] - K 2 F(r) = (1) 

in a spherical coordinate system to represent the field F and the position vector r can be traced back to 
as early as the 1860s, when Clebsch considered the scattering of an elastodynamic wave by a solid sphere, 
as recounted by Logan in 1965 [1]. In frequency-domain electromagnetics for homogeneous, isotropic, 
dielectric materials, this equation simplifies to the vector Hclmholtz equation 

V 2 F(r) + « 2 F(r) = (2) 

because of the constraint V • F(r) = 0. Credit for this reduction and application to the optical scattering 
response of a homogeneous, isotropic, dielectric sphere evidently goes to a series of papers by Ludvig 
Lorenz from the second half of the 19th century [1,2]. However, Lorenz's 1880 paper [3-5] has largely 
been superseded in popular memory by the 1908 paper of Gustav Mie [6], most probably because Lorenz 
used a pre-Maxwellian theory of light whereas Mie used the frequency-domain Maxwell equations and 
also applied his solution to the very practical problem of scattering of light by colloidal spheres. 

Solutions of Eqs. (1) and (2) in spherical coordinates are called vector spherical wavefunctions. For 
homogeneous, isotropic, dielectric materials, these wavefunctions involve n = ko^/e^-, where fcg is the 
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free-space wavenumber and e r is the (frequency-dependent) relative permittivity scalar. Extension to 
homogeneous, isotropic, dielectric-magnetic materials [7,8] is straightforward with the vector Helmholtz 
equation requiring only a simple redefinition: k = fco^/e r /i r , where /z r is the relative permeability scalar, 
also frequency dependent. Extension to homogeneous bi-isotropic materials, which require the solution 
of 

V x [V x F(r)] + ik (( - £)V x F(r) - k 2 Q (e r ^ r - C£)F(r) = (3) 

with £ and £ as frequency-dependent magnetoelectric scalars, came much later [9]; even so, the solutions 
of Eq. (3) are simple combinations of the solutions of Eq. (2) [9, 10]. 

Vector spherical wavefunctions for homogeneous anisotropic materials have not been found in closed 
form. A synthetic approach that requires finding vector cartesian wavefunctions [11] and subsequent 
transformation to spherical coordinates [12] docs not yield closed-form expressions for the vector spherical 
wavefunctions [13,14]. This is because a spherical coordinate system is generally incompatible with the 
cartesian nature of anisotropy. 

An exception is furnished by the Minkowski constitutive relations of a simply moving material that is 
an isotropic dielectric-magnetic material in its co-moving frame of reference. These constitutive relations 
can be stated as 

D(r) = e er[7l+(l-7)vv] • E(r) + /3v x H(r) (4) 
B(r)=/x /i r [ 7 |+(l_ 7 )vv] • H(r) - /3v x E(r) (5) 

where eo and are : respectively, the permittivity and the permeability of free space; e r , /U r , 7, and /3 
are real-valued scalars; the unit vector v is parallel to the velocity of the material; and 7 is the identity 
dyadic. By a combination of two transformations [15], the frequency-domain Maxwell equations in this 
material can be converted to frequency-domain Maxwell equations in a homogeneous, isotropic, dielectric- 
magnetic material [16]. The same algebraic manipulation is possible for a material whose constitutive 
relations are inspired by those of gravitationally affected vacuum and which can potentially be fabricated 
by properly dispersing electrically small bent-wire and other complex inclusions of different shapes and 
materials in some host materials [15]. 

We show in this paper that a similar combination of transformations can be fruitfully applied to the 
more general material described by the constitutive relations 

D(r) = e e r g ' A ' A ' g' 1 * E(r) - c 1 T x H(r) , (6) 
B(r) = /i 0j u r g • A ' A ' ST 1 ' H(r) + c^ 1 T x E(r) . (7) 

In these equations, cq = 1/^/eoMo is the speed of light in free space; e r ^ and \i r ^ are complex- 
valued scalars; S_ is either a rotation dyadic or a product of rotation dyadics; the diagonal dyadic A 
is positive definite with real- valued elements; and the gyrotropic vector T can have complex- valued 
components. For later use, we define the free-space impedance 770 = vMo/eo and the relative impedance 
Tj r {oS) = y/ /j, r (io)/e r (Lo). The dependences of the constitutive variables e r , /U r , A, and T on the angular 
frequency uj > have not been mentioned for notational compactness, and an exp(— iuit) dependence on 
time t is implicit. 

Accordingly, using a set of three transformations developed in Sec. 2, we determine closed-form 
expressions for vector spherical wavefunctions applicable to a material with the constitutive relations 
(6) and (7), as shown in Sec. 3. The derived vector spherical wavefunctions are used in Sec. 3.3 to 
formulate the T matrix for scattering by a three-dimensional object composed of the chosen material. 
Section 4 is devoted to the radiation field of an electric source current density distribution in the embedded 
material: a closed-form coordinate-free expression of the dyadic Green function is obtained, followed by 
the elucidation of its singularity behavior and its representation in a bilinear form involving the derived 
vector spherical wavefunctions. 
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2 Transformations 



Representing the position vector r = [x, y, z) = (r, 8, </>), and denoting the source current density by J(r), 
we aim to solve 

V x E(r) = icjB(r) , (8) 

VxH(r) = -iwD(r)+J(r), (9) 

when Eqs. (6) and (7) have been substituted therein. This is done after applying three transformations in 
succession so that we have to deal with effectively a homogeneous, isotropic, dielectric-magnetic material. 

2.1 First transformation 

Defining the auxiliary fields [17] 

e(r) = E(r) cxp (-ik T • r) , (10) 

h(r) = H(r) cxp (-ik T • r) , (11) 

and auxiliary source current density 

j(r)=J(r)exp(-iAor.r) , (12) 
and noting that ko = uj/cq, we transform Eqs. (8) and (9) to 

V x e(r) = iujfio^r g ' A ' A ' hi' 1 * h(r) , (13) 

V x h(r) = -iuje e r g ' A' A' ^ ' e ( r ) + -K r ) ■ ( 14 ) 

2.2 Second transformation 

Next, we define another set of auxiliary quantities 

e'(r')=^ 1 -e(r), (15) 

h'(r')=^ 1 -h(r), (16) 

j / (r')=r 1 *j( r )' ( 17 ) 
and coordinates r' = (x' , y', z') = (r', 9' , (/>') such that 

r'=^T (18) 

in order to transform Eqs. (13) and (14) to 

V' x e'(r') = iuj^ r A 'A ' h'(r'), (19) 

V' x h'(r') = -iu;e e r A ' A ' e'(r') + j'(r') , (20) 

because 

g- 1 • (Vx|) •§. = V'x|. (21) 

This transformation merely amounts to a rigid rotation of the coordinate system, and can be avoided by 
an appropriate choice of the (x, y, z) coordinate system. 



3 



2.3 Third transformation 



Finally, we make use of an affine transformation [18, 19] associated with the scaling of space as 

r = 4«r', (22) 

where f = (x, y, z) = (f, 9, </>), and define the final set of auxiliary quantities as 

e(r')=4-e'(f), (23) 

h(r') = A ' h'(r) , (24) 

j(r') = (adjA) .j'(f), (25) 

to obtain 



V' x e(r') = iu[iQ[i r ih(r') , 

V x h(r') = -iu>e e r le(r') + j(r') . 



(26) 
(27) 



where A is the determinant of A. Equations (26) and (27) are identical to analogous equations for a 
homogeneous, isotropic, dielectric-magnetic material. 



From these two equations, we get 

V' x [V' x e(r')] - k 2 e r ^ r A 2 e(r') = iw^oMr Aj(r') , 



V'x 



V' x h(r') 



k 2 e rf i r A 2 h(r') = V' x j(r') 



(28) 
(29) 



In the absence of a source current density, these equations are the same as Eq. (1). 



3 Field representation in a source-free region 



Homogeneous solutions of the Eqs. (28) and (29) in the spherical coordinate system were obtained by 
Clebsch [1]. Thus, in the absence of the source current density, 



6 ( r ')= E E EE [«l€J^+^MiM+ c iN« n (to') 

j — 1,3 s — e,o n—0 m= 



and 



oo n 



K r ') = E E E E [«&« L sL (kr') + /3<& n M« n (kr') + 7 «„ N«„ (fcr') 



j=l,3 s— e,o n—0 m— 

where the wavenumbcr 

k = k Q A^.yfJj^\ 

the standard vector spherical wave functions [20] 

1 



V [Y flmn (0',0')jn(fcr')] , 

v 



Y smn (6',cj>')h£Xkr') 

M«„ (fcr') = V x [r'y smn (0',0')in(fcr')] , 
M£> n (fcr') = V x [r' Y smn (9', 4>') h£\kr') 

N« n (fcr') = iv' x M«„ (fcr') ; 



(30) 

(31) 
(32) 

(33) 

(34) 
(35) 
(36) 

(37) 
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the tesseral harmonics 

Y emn (9', <j>') = P™ {cos 9') cos (rruf/) , (38) 
Y omn (6', 41) = P™ (cos 6') sin (m0') ; (39) 

jn( • ) is the spherical Bessel function of order n, hn ( • ) is the spherical Hankel function of the first kind 
and order n, P™( • ) is the associated Legendre function of order n and degree m; and the coefficients 
0') i.0') M JJ) «0) 

^smn, u smn: ^smn, f^-smri: t^sm 

variations of e(r') and h(r'). 



tj\ (~\ (j\ rj\ {j\ (j\ 
dsmn, bsmn, Csmn, «smn, Psmm and 7smn have to be determined from various stipulations on the spatial 



As the identities V x L^n (kr f ) = 0, V • M^L (fcr') = 0, and V • N^2n (kr') = follow from the 
definitions (33)-(37), substitution of Eqs. (30) and (31) in Eqs. (26) and (27) yields the simplifications 

<z«» = 0, aW n = 0. (40) 

Furthermore, as V x M^„ (fcr') = fcN^L (fcr') and V x N^„ (fcr') = fcM^„ (fcr'), the same substi- 
tution also yields the simplifications 

,>0) ,-j.O) 

B (j) = lc smn (j) = IQsmn , . . . 

r J srnn / \ ) I smn / \ • V / 

r] rj r (uj) 

Accordingly, Eqs. (30) and (31) can be rewritten as follows: 

oo n 

6 ( r ') = E E E E ^ M«„ (kr r ) + c«„ N«„ (fcr')] , (42) 

j = l,3 s—e,o n— 1 m— 

oo n 

^ = "^77^ E E E E [ C -» M -« ( fcr ') + (^')] ■ ( 4 3) 

Let us note here that the set of functions Mimn (fcr') and N^m™ (fcr') is complete [21]. 

3.1 Vector spherical wavefunctions 

Inverting the relations (23) and (24), we get 

oo n 

e'(r') =4-1. E E EE M &» * r ') + c &» (*=f ' ' r ')] . ( 44 ) 

j — 1,3 s— e,o n— 1 m— 

oo n 

h>') =~r^-4- 1 - E E EE [4iLM<iL (fc^ 1 . r') 

+€nNl(C 1,r ')] (45) 
from Eqs. (42) and (43). Next, an inversion involving Eqs. (15) and (16) yields 

oo n 

e(r) =g-A~ 1 ' E E EE [ 6 -« M«„ (fc^ 1 • | • r) 

j — 1,3 s—e,o n—1 m— 

+44 n NW„(fc4- 1 .|.r)l , (46) 



oo n 



i^NKC 1 -l-r)] , (47) 



r '° 7 ''' j=l,3s=e,on=lm=0 



from where we finally obtain 

E M = £ ' A' 1 ' E E E E M«„ (fc^ 1 • S • r) 



oo n 



j — 1,3 s— e,o n— 1 m— 

+4iL N«„ 1 • 5 • r)l exp (zfc r • r) , (48) 



oo n 



H(r) = --^-£-4- 1 ' E EEE^Mi^" 1 '!") 

'° ' r 3=1,3 s=e,o n=l m=0 

N[l (kAT 1 ' g ' r)] exp (*fc r • r) . (49) 

Thus, for the material described by Eqs. (6) and (7), the desired vector spherical wavefunctions in 
the spherical coordinates (r, 9, <j)) are 

m«„(r) = exp (zfc r • r) g • • M«„ (fc^ 1 • £ • r) (50) 

and 

n«„(r) = exp (*fc r • r) £ • A" 1 • N« n (fc^" 1 • £ • r) . (51) 
Equations (48) and (49) can then be compactly written as 



E ( r ) = E E E E [ b *rLn m&„ W + (r)] , (52) 

j — 1,3 s— e,o n— 1 m=0 

oo n 

H(r) =-— E E EE R™« to + to] " ( 53 ) 

/u " j=l,3s=e,on=ltn=0 

The completeness of the set of functions mim„ (r) and nimn (r) follows from the completeness of the set 
of functions (At') and (Ax') [21]. 

3.2 Closed-form expressions 

Without any loss of generality, we set 

: — t *) 

(54) 



g = L 

A = a~ x xx + oC x yy + zz 



where eta; > and a y > 0. In order to avoid extremely cumbersome expressions for the vector spherical 
wavefunctions uismn (r) and v^Xn ( r ), we also define wavefunctions A"t^„ (r) and A/"^ m (r) via 

m« n (r) = 4- 1 • M«„ W = R • -M«n (r) , (55) 

ni£„ (r) = ,r 1 * MiL W = I • AT«„ (r) , (56) 

with 

B = Ojxx + a y yy + zz . (57) 
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Using Eqs. (35)-(39), (50), and (51), we obtain 



M^ n (r) = exp [ik T ' r 



h{4>) 



sin 9 cos 6 Q smn {0, < 

f A (cf>)co S 2 8 + tf((b) S m 2 8 „ 

M6,<f>) Mmn{ 



(a x - oty) sin 9 sin <fi cos 72. 
i - (as - a y )cos(9sin0cos(/>7?. 







/ 2 (M) 



cos6>shi(/>cos0 Q smn (6>, 0) - h(<p)U smn (9,())) 



and 



Af«„(r)=exp [tt„r.r] 

"cos 2 6> + / 4 (<?(>) sin 2 6» 



x r 



+0 



kr 

K, ( j\kv) 
' hift>) 
Vi j) (fcr 

/Cl j) (fcr) 



fM - fK4>) 
MM) 



sin 6 cos 072, 



+ (a x — ay) sin sin </> cos </> Q 



smii V t 



/ 4 (0) cos 2 + /?(<£) sin 2 



/i(0) 

f ^(fcr 
fcr 

/i(0) 



/a(M) 



0) + (ctx — ctj/ ) cos sin <f> cos (/> <2 S 



sm W sm cos < 



/a 



■ cos 9 sin cos </> 1Z S 



h{4>) Qs 



Furthermore, 



R = [h{4>) sin 2 9 + cos 2 6>] rr + [/ 4 (c/>) cos 2 9 + sin 2 6>] 00 + [a x + a y - / 4 (</>)] 00 
+ [/ 4 (0) - 1] sin 6» cos 6» (r0 + 0r ) 

— (a^ — Qfy) sin <f> cos sin 9 (r(f> + 4>v J + cos 6> ^00 + 00^ . 
In the foregoing expressions, the following angular functions have been used: 



( al cos 2 i 



2 • 2 A 1/2 

a; sin 0) 



/l(0) = 

MM) = + [/ 1 2 (0)sin 2 6' + cos 2 0] 1/2 
/ 3 (0) = tan" 1 (^tan^ , 



U{<t>) 



a r cos 



The angle / 3 (0) must lie in the same quadrant as its argument. Additionally, 
jW(kr)=j n ikrf 2 (6,<f>)] , 

n+1 



/C«(fc 

Qsmn ( 



*r/a(M) 
= n(n + l)P 1 

cos 9 



jtfXkr)- J&ikr), j €{1,3} 



Mo,' 



/ 2 (M) 



h{4>)sm0 



V sro (0), 

u sm {4>) , 



i 



h{4>) sin < 



cos 



-(n + l)cos 6»P^ 



COS0 



Usm(4>) 



Vsm{,4>) 



M 

-sin [m/ 3 (0)] 
cos [m/ 3 (^)] 

cos [m/ 3 (0)] 
sin [m/ 3 (0)] 



(65) 
(66) 

(67) 



(69) 

(70) 
(71) 
(72) 



3.2.1 Special case: a x = a y — a 

When the diagonal dyadic A is simplified from biaxial (a x ^ a y ) to uniaxial (a x = a y = a), the functions 
fi to fi simplify considerably. The wavefunctions Al^ n (r) and A/"^„(r) then simplify as follows: 



M[ j i n {T) = cxp [iko r • r] (fcr) <j - f — -f sin 8 cos 9 Q 



a - I 



. a sin 2 + cos 2 i 



/b(0) 

A/"il(r) = cxp [ifc T • r 

fU) 



h{0) 

Qsmn 



x < r 



+6/ 



jj j) (fcr) q-1 
fcr / 2 (0) 



a- 1 



Sin COS 9 IZsmn (d, 



sm ft cos 



^ P sm „ (9, </>) + ICiP (fcr) M£ ti smn (9, , 

/5W 



^/C«(fcr) Q smn (9, 



(73) 



(74) 



Similarly, _B also simplifies to 

B = / 6 (0) rr + [a + 1 - / 6 (0)] 00 + a cf>4> + (a - 1) sin 9 cos 6» ( r6 + Or 



(75) 
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The following shorthand notation has been used: 

h{9)= + (a 2 sin 2 9 + cos 2 9) 1/2 , 
fe(6) = a sin 2 9 + cos 2 6 , 

Ji 3 )(fcr) = / l « [fcr/ 5 (0)] , 



Qsmn(@, 1 



^'W-J^W, j £{1,3}, 



n + 1 

= n(n + l)P, 

= < 
1 



COS 



/ 5 (0) 



cos (m^) 
sin (m<p) 

— sin (m<j 



a sin 9 \ cos(m0) 



asmtf 
(n + l)cos0P^ 



( n - m+ l)/ 5 (0)p™ 



cos 

W) 



cos (m<^) 
sin (m^>) 



(76) 
(77) 
(78) 
(79) 

(80) 
(81) 
(82) 



(83) 



3.2.2 Special case: a x = a y = 1 

Further simplification to A = J leads to B = 7, so that 



m« n (r) = M« n (r) = exp (*fc T • r) M« n (fcr) , 
n«„(r) = ATi£„(r) = exp (*fc T • r) N«„(fcr) . 



(84) 
(85) 



3.3 T matrix for scattering by a 3D anisotropic object 

For an example of the use of the derived vector spherical wavcf unctions, let us now consider scattering 
by a 3D object that occupies the region V and has a surface S. The material occupying V obeys the 
constitutive relations (6) and (7). The coordinate system is chosen so that its origin lies inside V and 
the conditions (54) hold. The external region is vacuous. 

Suppose that an electromagnetic field is incident on the scatterer. This incident field is purely arbi- 
trary except that its source must lie outside a minimum sphere centered at the origin and circumscribing 
V. It can then be represented as [22] 



oo n 



E mc (r)= J2 EE D ™ €„Mil(fcor) + 5(lN(l(fc r) 



s—e,o n—1 rn—0 

oo n 



s—e.o n—1 m—0 



H mc (r) = Yl EE °™ ^iNlK^+BlMl^r) 

m ^ ^ ^— ' 



(86) 
(87) 



where the coefficients A^] ln and B^rL n are presumed to be known and the normalization constant 



(2 - 5mo) 



2n+l (n-m)\ 



An(n + 1) (n + m)\ 



involves the Kronecker delta S pq . Outside the circumscribing sphere, the scattered electromagnetic field 
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is represented by the sums 

E sca (r)= EE D ™ 4 3 2 n M( 3 2jfcor) + S(lN( 3 „)„(fc r) 



OO 71 



s— e,o ri— 1 m=0 

oc n 



H sca (r) = -- E EE D ™« ^inN^JfcoO+S^LM^Jfcor) 



(89) 
(90) 



s— e,o n— 1 m— 



where the unknown coefficients AUn and B S m n have to be determined. The underlying presumption 
here is that the source of the incident field lies so far away from the scattering object that any action of 
the scattered field on that source is vanishingly small. 

Application of the Ewald-Oseen extinction theorem followed by some analytic-continuation argu- 
ments leads to the following four relationships [22,23]: 

A&n = - V / ^ {^^^ X Emt(rp)] ' N — (k ° rp) 



D 



(i) 



ir]„ [n(r p ) x H int (r p )] • M^„(fc r p )} , 

= J s d \ {[ n ( r p) x Ei„t(r p )] • M^„(fc r p ) 

+ir]a [n(r p ) x H int (r p )] • N^> n (/c r p )} , 

4 3 in = ^ J/*v X Emt(rp)] ' N — ( fcor ^ 

+%o [n(r p ) x H int (r p )] • M^Jferp)) , 

= v i d2rp { [n(rp) x Emt(rp )] ' M ^"(*° r p) 

[n(r p ) x H mt (r p )] • N« n (fc r p ) j . 



(91) 



(92) 



(93) 



(94) 



Here, n(r p ) is the unit outward normal to S at r p <G S, whereas Ei„t(r) and Hi„t(r) are the fields excited 
inside the scatterer. 



A representation of the internal fields emerges from Eqs. (52), (53), (55), and (56) as 



oo n 



^ = 1 * E E E [ b $L M *™ (r) + cV n JtfVL (r) 

e,o n—1 m— 

oo n 

-i • E E E [ c -« • M -« ( r ) + €„M2 n (r) 



s — e,o n—1 m— 

because the internal fields must be regular at the origin OeV. 



(95) 
(96) 



Substitution of Eqs. (95) and (96) in Eqs. (91) and (92) leads to the matrix equation written com- 
pactly as 

, , , 

srnn.s'm'n' 



r a (i) ' 




o(l) 
J-'smn 









K 



smn.s' in' rv 



J (1) 

smn,s'm' n' 




r , , i 

s'm'n' 


, , , 

smn,s m n 




<£> , , 

L s m n J 



y(i) 



6^ , , 

s m n 



(97) 



10 



likewise, Eqs. (93)— (96) yield written compactly as 



r a {3) ' 




r / (3) 

smn^s' m'n' 


r(3) 

smn.s'm'ji' 




r 6 ( P , , i 

s m n 


B (3) 




K {3) , , , 


r(3) 
smn..s'rn'n' 




c ( ? , , 

L s m n J 



y(3) 



, , 

s'm'n' 



s' nvn' 



These relations require the computation of the integrals 

,2 



I 



GO 



sni it .s' rn' n 



J, 



go 



sin ii. s' Hi/ a 



K 



GO 



sin 1 1 ,s m n 



- GO 



smn.s' m' n 



jf d 2 r p (Nf)„(t,r p ) • {n(r p ) x [g • M« , n ,(r p )] } 
ML(¥ P ) • {n(r p ) x [l'A/^, n ,(r p )]}) , 

| s d 2 r p (NW m (fc r p ). {n(r p )x • « , n , (r p )] } 
iMW„(fcor P ) • |n(r p ) x [fl • M&, n , (r p )l V) , 

,'m'n' = - V l ^ ( M — ( fcor ^ * X [l * *tfL»'("p)] } 

^NKfcor,) • {n(r p ) x • M^ m , n , (r p )] }) , 

_!M / ^ ( M W n (fc r p ) • {n(r p ) x • A#>>, n ,(r p )] } 



t ( fc or p ) • |n(r p ) 



'm'n'( r p) } 



where j G {1, 3} and I = j + 2 (mod 4) G {3, 1}. 



(98) 



(99) 



(100) 



(101) 



(102) 



Equations (97) and (98) yield a matrix relationship between the incident-field coefficients and the 
scattered-field coefficients that can be expressed as 



(103) 



r ^ (3) " 




' A (1) 1 




= [T] 




o(3) 




o(l) 
^ sinn 



where the T matrix 

[Tl = - \y {3) ] _1 

completely characterizes the scattering response of the 3D anisotropic object. 



(104) 



The calculation of the T matrix of a 3D object made of a material obeying the constitutive relations 
(6) and (7) along with (54) is a major exercise that we intend to take up shortly. Let us, however, mention 
that our derivation procedure can be extended to cover the possibility that the external region is occupied 
by a homogeneous isotropic material which can be dielectric- magnetic, chiral, or even biisotropic [24] . 
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4 Radiation fields 

Returning to Eqs. (8) and (9), we are now interested in finding the electric field radiated by a source 
current density distribution. Accordingly, 

E(r) = iuifj,QH r J J J <Py p Q{y,y p ) • J(r p ) , (105) 

where G(r, r p ) is the dyadic Green function with r as the field point and r p as the source point. 

4.1 Dyadic Green function: Coordinate-free form 

Our starting point is the particular solution of Eq. (28), which is known to be [25] 

e(r')=i^ of x r A [ j [ d 3 r' p ~g(r'y p ) >](r' p ), (106) 



where 



Inverting the relations (23) and (25), we get 

e'(t>) J J |d 3 r;|'(r',r;).j'(r;) (108) 

from Eq. (106), where 

£'(r',r p ) = A 3 A-' '-giA- 1 -r'^A- 1 'r^'A- 1 (109) 

i n , , i exp \ik\A~ 1 • (r' - r' ) |1 

= A 3 A • (/ + k~ 2 A • V V • A) • A — — — ■= ^ iLll (no) 

= v = = =' = A-k\A~ 1 ' (r' -r' p ) | 

= A 3 (A + fc V V ) 1 = ^ P-^-, 111 



and the identity A T = A nas been used. 



Next, an inversion involving Eqs. (15) and (17) transforms Eq. (108) to 

e(r) =iujn n r J J J d 3 r p g(r,r p ) • j(r p ), (112) 

where 

g(r,r P ) = g'^(g'r,g'r p ) • g' 1 (113) 

„ „ „ , exp [ifclA -1 • S • (r - r„) |1 

= A 3 S • (A + k S • VV • S ) • S 1 '= ^— 114 

= V = = = ' = AirlAT 1 ' g ' (r - r p ) | V ; 

= A 3 (S • A~ 2 • S-i + k~ 2 S 2 • VV • S->) eXP ^ 1 ^- (r "^ )l] . (115) 
V = = = = = ' ^lAT 1 • g • (r - r p ) | 1 j 

Finally, the use of Eqs. (10) and (12) in Eq. (112) yields Eq. (105) with 

£(r,r p ) = exp[ifc r • (r-r p )]g(r, r p ) (116) 

= ,4 3 exp [ifc r • (r - r p )] (S ' AT 2 ' gT 1 

9 9 9% ex P Rfcl^l^ 1 * S • (r - r„) 11 

+k~ 2 S 2 • VV • S- 2 ) 1 = ^— 51! H7 

= 7 47TI4- 1 (r-r p )| 1 ; 
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in a coordinate-free form. 



4.2 Dyadic Green function: Singularity 

When the source point r p and the field point r coincide, the dyadic Green function exhibits a singularity. 
In order to estimate this singularity, the source point is made the centroid of an electrically small convex 
region bounded by the surface S^, which is then shrunk in volume isotropically [26]. This process is 
used to determine the depolarization dyadic D defined from G(r, r p ) via 

lim| J jf d\g(r p ,r q ) • b(r,) ~ £ • b(r p ) , (118) 

where the surface 

S e = {r q \r q {n i ) = r p +lg e -u i (n i ), 0^ e [0,4tt]} . (119) 

Here, the positive scalar £ is a linear measure of the volume of ; the dyadic ^ delineating the shape 
of has a determinant equal to unity, while ii£ and fli denote, respectively, the radial unit vector and 
the solid angle in a spherical coordinate system with its origin at r p 6 V|. The depolarization dyadic D_ 
depends on the choice of ^ [26] . 

In the present instance, regardless of the choice of two different approaches [27] to determining 
D_ arc possible as follows: 

(a) Determine the analogous depolarization dyadic d arising from g(r' ,r' p ), and then reverse the three 
transformations described in Sees. 2.1-2.3 to obtain D_. Implementing this indirect approach, we 
find that Eq. (107) leads to 

d=-k~ 2 L, (120) 



whore 

XlqYq UqY T j 

'St ! 47r|r ? -r p |- 

and u q is the unit outward normal to at r q £ S^. The dyadic L is symmetric [26]. If is either 
spherical or cubical, L = (1/3)1. For ellipsoidal V^, the dyadic L can be cast in terms of elliptic 
integrals that have to be evaluated numerically [28,29]. Using Eqs. (109) and (120), we get the 

i = -k- 2 A 3 A^ ' L ' AT 1 (122) 

arising from g'(r',r' p ), and finally 

£> = -k' 2 A 3 g -AT 1 'L 'AT 1 '§T X (123) 
emerges from the use of Eqs. (113) and (116). 

(b) Directly obtain D_ for the material described by the constitutive relations (6) and (7), using the 
method of Michel and Weiglhofer [30] . This direct approach was formulated for ellipsoidal and 
yields [30] 



A 2 [ 2 * r UZ 1 ' pp ' C/7 1 
D = -— j / d<p p d& p wa.tio r-^ =* , , (124) 

where p = (x cos </? p + y sin <^ p ) sin $ p + z sin z? p . 
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For spherical V £ , the right side of Eq. (123) simplifies to -(A 3 /3k 2 )S_ 'A 2 • S, but the integral 
on the right side Eq. (124) has still to be expressed in terms of elliptic integrals [31]. Therefore, the 
depolarization dyadics (123) and (124) — obtained by two different approaches — are clearly distinct in 
general. The difference arises because, for the two approaches, 

(i) \/i was chosen in different spaces: r' in the indirect approach and r in the direct approach; and 

(ii) the singularity was considered for different fields e in the indirect approach and e (or E) in the 
direct approach. 

Knowledge of the depolarization dyadic will be useful in setting up the Maxwell Garnett and the Brugge- 
man formalisms to predict the effective constitutive parameters of a particulate composite material 
containing electrically small particles made of the chosen material [32] . 



4.3 Dyadic Green function: Bilinear expansion 



Computation of the field radiated by a source can be greatly convenienced by the formulation of a 
blinear expression of the dyadic Green function G(r,r p ), because the bilinear expression is a sum over 
the products of the derived vector spherical wavefunctions of the source point and the field point [33,34]. 
It is applicable only if the source point and the field point are situated at different distances from the 
origin. 

We begin with the known bilinear expansion [22,23] 

■ L. oo n 

~g(r',r' p ) = ^EEE fl ™ [MlKlMlKKl + NiKJNlK)] , (125) 



s— e,o n— 1 m— 



where 



r' > r' p 
r' < r' 



By following the process described in Sec. 4.1, we get 



r> >r' p 
r' < r' 



ikA 3 



oo n 



, 'A'' • E E E D ™ [MifnnikA- 1 ' rgMWjfc^" 1 • r'<) 

s— e,o n— 1 m— 

+N[i(^ 1 ' r^UkA- 1 • r'J 



A- 1 



and 



S0,r p ) 



ikA 3 



exp [-ik T ' (r> + r<)] 



OO n 



E EE D <™ m^„(r > )mW„(r < ) + n^ n (r > )nW„(r < ) 



s— e,o n— 1 m— 



where 



r> 



r > r„ 



r > r p 
r < r v 



(126) 



(127) 



(128) 



(129) 



Finally, from Eqs. (116) and (128), we get the desired bilinear expansion 

ikA 3 



G(r,r p ) 



7T 

oo n 



exp [-ik T • (r> + r< - r + r p )] 



E EE °™ m( 3 J l Jr > )m(^(r < ) + n( 3 J l Jr > )nW m (r < ; 



s—e,o m— 



(130) 
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Application of the conditions (54) simplifies the foregoing expression as follows: 
G(r, r p ) = exp [-ik T • (r> + r< - r + r p )] 

— 7T 

oo n 
s— e,o n— 1 m— 

Equations (130) and (131) will be helpful for analyzing aperture antennas [33,34] and radomes made of 
the chosen material [35,36]. 



5 Concluding remarks 

The derivation, in closed-form, of vector spherical wavefunctions for an orthorhombic dielectric-magnetic 
material with gyrotropic-like magnetoelectric properties, enables the electromagnetic scattering response 
of an object made of such a material to be comprehensively and conveniently characterized. In partic- 
ular, these wavefunctions provide the building blocks for (i) the T-matrix formulation which yields a 
description of the field scattered by a 3D object composed of the chosen material, and (ii) the dyadic 
Green function which yields a description of the radiation field of a source current density distribution 
embedded in the chosen material. 

Acknowledgments. AL thanks Craig F. Bohren for a discussion regarding the antecedents of the 1908 
paper of Mie [6] , and the Charles Godfrey Binder Endowment at Penn State for ongoing support of his 
research activities. 



A4i 3 2jr > )A4Wjr < )+AAWjr > )AAW, l (r < ) > B. (131) 
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